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Executive  Summary 


This  report  documents  the  results  of  two  tasks.  The  first  task  was  the  continued  study  of  the  Raman 
dynamic  secular  equation  of  cubic  crystals,  which  relates  the  Raman  frequency  shifts  and  the  surface 
stress  field.  The  Raman  frequency  shifts  under  a  general  two-dimensional  stress  field  were  formu¬ 
lated.  The  second  task  was  an  experimental  effort  using  Raman  spectroscopy  to  measure  the  fre¬ 
quency  shifts  in  a  single-crystal  structure  that  can  be  translated  into  surface  stress  fields. 

In  the  first  task,  the  three  Raman  frequency  shifts  under  bi-axial  stresses,  x„,  tjj,  and  a  shear  stress 
in  a  local  coordinate  system  for  cubic  crystal  systems,  were  determined  analytically.  The  relation¬ 
ship  between  the  three  Raman  frequency-shift  values  and  the  three  stress  components  is  unique. 
Therefore,  for  a  set  of  Raman  frequency  shifts,  the  three  stress  components  will  be  uniquely 
determined. 

In  the  second  task,  Ramon  frequency  shifts  were  measured  using  single-crystal  silicon  strips  under 
uniaxial  stress  along  either  the  (100)  or  the  (110)  orientation  of  the  crystal  system  in  a  test  fixture  to 
test  the  derived  expressions.  The  measured  frequency  shifts  based  on  the  calculated  stresses  agree 
well  with  the  theory  derived  from  the  first  task,  thus  validating  the  model. 


V 


1.  Objective  and  Approaches 


1 .1  Objective 

The  objective  of  the  work  is  to  develop  innovative  experimental  techniques  for  determining 
mechanical  properties  for  MEMS  elements.  Specific  tasks  are: 

(1)  Derivation  of  relationships  between  Raman  Shifts  and  applied  two-dimensional 
surface  stresses,  and 

(2)  Development  of  innovative  experimental  techniques  to  determine  mechanical 
properties  important  to  MEMS  components. 


1 .2  Approaches 

For  the  duration  of  this  program,  two  tasks  were  performed  to  achieve  the  objective.  The  tasks 
included  an  analytical  evaluation  of  the  Raman  dynamic  frequency  equation  and  an  experimental 
testing  method  to  detect  stresses  under  external  loads  applied  primarily  to  single-crystal  silicon  strip 
specimens.  The  Raman  dynamic  frequency  equation  was  examined  in  more  detail,  and  the  frequency 
shifts  under  general  two-dimensional  stress  states  were  determined  analytically.  Because  of  the 
capability  of  the  Raman  technique  in  addressing  very  localized  stress  (spatial  resolution  1-2  |xm), 
potential  stress  non-uniformity  can  be  identified.  An  observed  stress  non-uniformity  in  an  otherwise 
theoretically  uniform  stress  region  will  lead  to  a  conclusion  of  material  inhomogeneity.  In  turn,  this 
understanding  should  lead  to  improved  designs  of  MEMS  devices  using  validated  solid  mechanics 
principles. 
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2.  Introduction 


The  use  of  Raman  spectroscopy  is  becoming  known  as  a  diagnostic  in  mapping  stress  in  composites, 
semiconductors,  and  ceramics  over  the  past  few  years  and  needs  to  be  refined  for  space  applications. 
For  example,  Raman  has  been  used  with  some  success  to  study  relative  stress  at  interconnects,  near 
local  oxidation  structures,  vias,  and  gate  oxide  structures,  and  to  study  adhesion  effects  between  thin- 
film  layers  as  affected  by  stress.  The  technique  is  also  applicable  to  stress  migration  (SM)  and 
electro-migration  (EM)  (two  phenomena  observed  during  life  and  reliability  testing  of  interconnects) 
and  therefore  to  circuit  reliability. 

Raman  spectroscopy  is  a  non-contact  diagnostic  technique.  The  Raman  effect  is  observed  as  either 
the  gain  or  loss  of  photon  energy  as  a  result  of  photon/phonon  scattering  from  molecules  or  crystal 
lattices.  The  observed  transitions  are  an  indication  of  resonance  in  the  vibrational  energy  of  a  mate¬ 
rial.  As  the  material  is  strained,  the  energy  state  in  the  material  is  altered  and  thus  the  frequency  is 
changed.  The  spatial  resolution  of  Raman  measurements  can  be  made  very  high,  on  the  order  of  1 
|xm.  It  is  very  useful  in  detecting  the  surface  stresses  in  either  single-crystal  or  polycrystalline  sys¬ 
tems.  The  Raman  spectroscopy  technique  is  used  here  to  demonstrate  the  detection  of  the  residual 
stresses  in  a  frequently  used  MEMS  cantilever  beam  structure  and  in  single-crystal  test  articles. 

In  this  report,  the  Raman  secular  equation,  which  determines  the  Raman  frequency  shifts  under  an 
applied  stress  field,  is  first  discussed.  The  experimental  results  on  the  Raman  frequency  shifts  using 
various  types  of  specimens  and  different  stress  fields  are  also  discussed. 


3.  Raman  Spectroscopy 


3.1  Raman  Secular  Equation 

The  dynamic  equations  for  the  optical  modes  in  diamond-like  cubic  crystals  have  the  form'’^ 


mu  =  H  Kf%i 


+ 


Z 

klm 


ik 


'nim^k 


(1) 


where  Uj  is  the  i*  component  of  the  relative  displacement  in  the  unit  cell;  m  is  the  reduced  mass  of  the 
atoms;  K./’  =  m  (o/  is  the  effective  spring  constant  in  absence  of  strain;  (3K„/9ti,Jt1|^  =  K^Tii^  is  the 
change  in  the  spring  constant  due  to  an  applied  strain  Ti,^;  and  i,  k,  1,  and  m  refer  to  x,  y,  or  z  in  an 
xyz  coordinate  system. 

are  elements  of  a  fourth-order  tensor.  They  are  the  second  derivatives  of  the  potential  energy 
with  respect  to  the  displacement  fields  u^.  Their  magnitudes  depend  on  the  types  of  crystals.  The 
crystal  symmetry  rules  also  apply  to  them.  For  a  cubic  crystal,  for  example,  there  are  only  three 
independent  components,  “  ^^2222  ^^333  ^  P,  ^1122  ^1133  ^^2233  ^  ^tlclK^2^2  ^1313  ^^323 

=m  r.  Consequently,  one  obtains  the  following  secular  equation,  whose  solutions  yield  the 
frequencies  of  the  optical  phonons  in  the  presence  of  strain.  The  constants  p,  q,  and  r  are  measured 
values  for  each  type  of  crystal. 


[p^xx  “*'^(^yy '*'^zz) 

2rexy 

2rexy 

Ip^yy  ^^(^xx  ■*‘®zz)~^J 

2r£vz 

r  /  ^  \  1 

=0 

(2) 

2r8xz 

2r£  y2 

IP^zz  “^^l^xx 

where  the  are  the  strain  components,  and  X  =  £2^-co/  represents  the  frequency  shift.  Since  X  is  so 
much  smaller  than  the  shifted  frequency  Q.  can  be  approximated  by  £2  =  +  Xy2co^and  A£2  =  £2 

-  (0,  =  A/2(0„. 

It  is  seen  from  Eq.  (2)  that  under  an  incident  monochromatic  light  or  laser,  there  will  be  up  to  three 
Raman  frequency  shifts.  The  shifts  are  a  function  of  the  elastic  strains  that  are  uniquely  determined 
from  the  applied  stresses.  The  shifts  associated  with  some  special  stress  states  are  discussed  in  the 
next  section. 

3.2  Frequency  Shifts  Associated  With  Two-Dimensional 
Stress  Field  and  Examples 

The  strain  components  e,j  are  related  to  the  stress  components  i:~  by  the  following  relationship: 
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{£}  =  [S]{X} 


(3) 


where  [S]  is  the  compliance  matrix.  Its  elements  are  components  of  a  fourth-order  tensor.  For  a 
material  that  is  cubic  such  as  silicon,  the  [S]  has  the  following  form: 


^1 

®21 


0 

0 


S^2  ^3  0  0  0 

S22  ^23  0  0  0 

%2  ^3  °  °  ° 

0  0  S44  0  0 

®  °  *3  ^5  0 

0  0  0  0  Sgg 


(4) 


The  Raman  secular  equation  reduces  to  a  cubic  equation  for  X  for  given  values  of  elastic  strains  asso¬ 
ciated  with  the  corresponding  stress  field.  The  solutions  are  designated  as  A.,,  X^,  and  A3. 

The  Raman  frequency  shifts  associated  with  three  different  stress  fields  were  analyzed;  case  1  is  a 
uniaxial  stress  field  in  the  x-direction;  case  2  is  an  equal  biaxial  stress  field  in  the  x-  and  y- 
directions;  and  case  3  is  an  equal  tension  and  compression  in  the  x  and  y-directions.  Case  3  also 
corresponds  to  a  pure  shear  stress  field.  Since  Raman  is  a  surface  observation,  the  frequency  shifts 
are  only  related  to  three  stress  components  at  the  observed  surface,  one  shear  stress  and  two  normal 
stress  components.  After  solving  Eq.  (2),  the  unique  frequency  shifts  associated  with  each  case  are 
shown  in  Table  1. 


The  Eigenvalues  from  Eq.  (2)  can  be  treated  as  a  linear  combination  of  the  Eigenvalues  of  two 
decomposed  stress  fields  for  any  given  stress  field,  one  associated  with  the  hydrostatic  stress  field 
and  the  other  associated  with  the  deviatoric  stress  field.  The  hydrostatic  and  deviatoric  stress  fields 
are  defined,  respectively,  as: 


Nh  =  'to  =  ('C„+t^+T33)/3 

=  (5) 


where  5.  is  the  Knoecker’s  delta. 

Using  and  as  the  stress  fields,  the  corresponding  strain  fields  are  then  computed.  The 
three  Eigenvalues  for  {Xjj}^are  identical.  The  three  Eigenvalues  for  {Xy}^may  be  distinct.  Let  A^  be 

Table  1.  Raman  Frequency  Shifts  Corresponding  to  Various  Liniiting  Stress  Cases 


Frequency  Shifts 


Case  1 

A£l,=  i:.(pS„+2qS,J/2<D. 

AQj=  TjpS„+q(S„+S„)]/2(0. 

Case  2 

AQ,=An,=  Tjp(S„+SJ 

An,=  T„[pS„+q(S„+S„)]/(0„ 

+q(S„+3SJ]/2(D„ 

Cases 

Xo(P-q)(S„-S,j)/2<B, 

0 

il 

AQ,=  -T„(p-q)(S„-S„)/2to„ 
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the  Eigenvalue  associated  with  and  let  and  X^  be  the  three  Eigenvalues  associated  with 

{%}.■ 

By  observation,  we  have 

\  =(  X^+Xj  +  XgVS 
X„=(X,-A,)/3+(X,-X3)/3 
K=iK-\V^HK-KV3 
X,,=iX,-X,)/3+{X,-X,)/3 


or 


x,=x^+x^, 

X2=X^+X^ 

The  frequency-shift  values  in  terms  of  hydrostatic  and  deviatoric  components  are  shown  in  Table  2. 
The  frequency  shifts  are  independent  of  the  frequency  of  the  incident  light.  It  is  seen  that  for  an 
equal  biaxial  stress  field,  the  frequency  shifts  under  hydrostatic  stress  components  are  twice  those  of 
the  uniaxial  stress  field.  For  the  equal  tension-compression  case,  there  is  no  shift  associated  with  the 
hydrostatic  stress  because  it  is  zero.  There  are,  however,  two  deviatoric  stress-related  shifts  on  either 
side  of  the  reference  frequency  with  the  same  magnitudes  in  frequency  difference.  The  frequency 
shifts  corresponding  to  a  general  stress  field  can  be  obtained  numerically. 


Table  2.  Raman  Frequency  Shifts  Corresponding  to  Hydrostatic  and  Deviatoric 

Components 


Frequency  shifts  under  hydrostatic 
stress  components 

Frequency  shifts  under  deviatoric  stress 
components 

Case  1 

^h=To{P+2q){S„+2S,2)/6(Oo 

Afl„=-Co(P-q)(S„-S,j)/3(0„ 

AQ^=A£l*,=  -An„/2 

Case  2 

Afl„=T.(P+2q)(S„+2SJ/3«), 

=  x„(p-q)(S„-S  J/6co„ 
An«=-2An„ 

Case  3 

^d,=t«(p-q)(s„-sj/2(Oo 

Aii«=-to(p-q)(s„-s,j/2(o„ 

An^  =  o 

3.3  Frequency  Shifts  Associated  With  Two-Dimensionai  Stress  Fieid 

The  two-dimensional  stress  field  is  appropriate  to  a  surface  where  stress  perpendicular  to  the  surface 
is  absent.  This  assumption  applies  to  most  surfaces.  Consider  this  general  two-dimensional  stress 
field:  x„,  x^j,  and  x,j  where  x„  and  x^j  are  normal  stresses,  and  x,j  is  the  shear  stress  with  respect  to  a 
Cartesian  coordinate  system  with  axes  1, 2,  and  3.  Further,  let  the  1, 2,  and  3  axes  coincide  with  the 
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principal  axes  of  the  cubic  crystal  material.  If  axis  3  does  not  coincide  with  the  principal  axis  of  the 
crystal,  there  will  be  six  stress  components  present.  This  makes  the  analysis  much  more 
complicated. 

Let  the  three  stress  components  be  related  by  the  parameters  a  and  p  as  follows; 

t„  =  T,  =  andT,3  =  pT  (6) 

Eqs.  (2),  (4),  and  (5)  combined  with  =  Q  -  co^  =  XI2L(Si^  give  the  following  three  frequency  shifts: 

AI2,  =  (T/4®J{p(l+a)(S„  +  S,,)  +  q(l+a)  (S„  +  3S,,)  +  [(l-a)\S„  -  S.,)Vq)' 

+  4(rSj)^]''^}, 

(t/4®^{p(l+a)(S„  +  S.,)  +  q(l+a)  (S„  +  38,3)  -  [(l-a/(S„  -  Sj\p-qf 
+4(rS,P)1'"}, 

Ml,  =  ('t/2®^{p(  l+a)S,,  +  q(l+a)  (S„  +  S,,) }  (7) 

The  hydrostatic  and  deviatoric  components  of  the  frequency  shifts  can  be  expressed  as: 

AQ,  =  (t/6®„){(p  +  2q)(l+a)(S„  +  2S,,)} 

AH,.  =  (x/12®J{(p  -  q)(l+a)(S„  -  S.,)  +  3[(p-q)^  (l-a)^(S„  -  S,/+  4(rS^P)^]''^} 

AO,,  =  (T/12mJ{(p  -  q)(l+a)(S,.  -  S„)  -  3[(p-q)^  (l-a)\S„  -  S„)V  4(rS,,P)^]“} 

A0,3  =  -(t/6®,){(p  -  q)(l+a)(S„  -  S,,)}  (8) 

The  Eigenvectors  associated  with  the  three  Eigenvalues  and  X,  are 


1 1 

-  A2 

"o' 

1 - 

0  > 

1 _ 

1 

V2  = 

1 

0 

1 

V3  = 

0 

_1_ 

Vl  +  A?  ’ 

<1 

+ 

(9) 


where 


^  _(l-cx)(Sii  -Si2)(p-q)/2-D  ^  _ _ rS44P _ 

^  rS44p  ’  2  (i_a)(Sii-Si2)(p-q)/2  +  D  ’ 

D=l^(l-a)\s„-Si2)^(p-q)2+4(rS44P)2 
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By  observation,  we  note  that  A,  =  -A^  and  all  three  Eigenvectors  are  mutually  orthogonal. 


3.4  Selection  Rule 

The  three  frequency  shifts  associated  with  the  stress  field  are  not  necessarily  observable  simultane¬ 
ously  at  a  specific  orientation.  Their  detectability  depends  on  the  combination  of  the  polarization  of 
the  incident  light  and  scattered  light.  For  cubic  crystals  such  as  silicon,  the  Raman  tensors  in  the  1, 
2,  and  3  directions  are 

'O  0  0]  [0  0  d' 

Rl=  0  0  d  ,  R2=  0  0  0, 

0  d  oj  [d  0  0_ 

The  Raman  tensors  corresponding  to  three  Eigenvectors  are 


'0  d  o' 

R3=  d  0  0  (11) 

0  0  0_ 

where  d  is  the  lattice  constant.  Its  magnitude  is  of  the  order  of  10  nm. 

Let  Cj  and  be  the  incident  and  scattered  light  with  polarization  of 

'ii 

ei=lk^ 

where  are  directional  cosines,  and  (//j  J  +  (//2  f  -i-  (//3  f  =  1 ,  k  =  i,s,  and  j  labels  the  frequency 
shifts,  A£2j  in  Eq.  (7). 

The  scattering  efficiency,  I,  is  defined  as  CL  le*.  R.  e'/  where  C  is  a  constant.  Accordingly, 

11  =  C  [ijsUh  Ai  +4 

12  =  C  [/,3  (I si  +  ls2^2  +  (13) 

13  =  C[{iyi+lll^2)df 


■  0  0  Ajd]  [0  0  d  ' 

Rl=^=i=  0  0  d  ,  R2=-r=  0  0  A2d  , 

Aid  d  0  J  [d  A2d  0  _ 


0  d  0 

R3=  d  0  0  (10) 

0  0  0_ 
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It  is  seen  from  Eq.  (13)  that  for  the  cases  that  correspond  to  backscattered  light  from  the  <001>  sur¬ 
face,  only  z-polarized  phonons  can  be  observed.  In  addition,  when  the  polaiization_^gles  of  inci¬ 
dent  light  and  the  scattered  light  are  complements  to  each  other,  the  expression  e^^  R3  e/  is 
maximized. 


3.5  Determination  of  Stress  Field  with  Known  Raman  Frequency  Shifts 

The  existing  literature  on  stress  measurements  using  Raman  spectroscopy  has  focused  on  Raman 
measurements  or  observation  when  the  applied  stresses  are  knovm  for  calibration  purposes.  Our 
challenge  is  the  determination  of  a  unique  stress  field  t„,  and  x,,  once  the  Raman  frequency  shifts 
are  measured. 

By  examining  Eq.  (7)  or  (8),  one  finds  that  the  three  frequency-shift  expressions  are  not  linearly 
independent  and  hence  cannot  be  used  to  uniquely  determine  the  stress  field.  Therefore,  one  more 
equation  is  needed.  In  backscatter  mode,  one  can  select  the  incident  and  scattered  polarization  as 


o' 

'o' 

T 

0 

,  eg  — 

1 

and  65  = 

0 

_0_ 

_0_ 

Inserting  these  polarizations  into  Eq.  (13),  one  obtains 

Il  =  Cd^/(l  +  Af), 

l2  =  C(A2d)^/(l  +  A^) 

using  the  first  e^,  and 

(15) 

Il  =  C(Ai  d)^/(l-t-A]), 

l2  =  Cd^/(l  +  A|) 

using  the  second  e^. 

(16) 

By  taking  the  ratio  of  I,  and  I^  in  Eq.  (15)  or  (16),  one  obtains  A,  or  A^.  a,  P,  and  x  can  then  be  deter¬ 
mined  from  Eq.  (7)  as 


X  =  (i)o 


AQ3  (Afl^-AQ^Xl-A?) 

pSi2+^(Sn+Si2)  (5ij-5j2)(p-9)(l+Af) 


a = 1 

(Sii-Si2)(p-9)(1+Af)x 


(17) 


p_  -Si2)(p-q)Ai 

(l-A?)rS44 
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It  is  noted  that  the  same  information  can  be  obtained  by  using  either  the  first  e^  or  the  second  e^  with 
Aj  =  -A,.  Also,  The  following  limits  of  A,  and  A^  exist: 

Aj ,  A2  e  (-1, 1),  Pe  (-00, 00) 

Asp->-~,  Ai-»-l,  A2^+1 

Asp->  Ai-^+1,  A2->-1 

At  p  =  0,  Aj  =  A2  =  0 

A  plot  of  A^  versus  P  for  values  of  a  =  -1,  -0.5, 0, 0.5,  and  1.0  is  depicted  in  Figure  1.  Also,  the 
analysis  has  shown  that  for  ipi  greater  than  1,  the  sensitivity  becomes  poor  for  a  given  value  of  a. 


Figure  1.  Plot  of  P  versus  Aj  for  single-crystal  silicon  (a=  -1,  -0.5, 
0, 0.5,  0.95). 
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4.  Raman  Stress  Experimental  Setup  and  Measurements 


4.1  Specimen  Description 

Experimental  test  calibration  specimens  were  made  of  single-crystal  silicon.  The  specimen  was  in  a 
strip  configuration.  The  strip  specimen  was  clamped  into  the  holder  at  its  two  end  regions.  The 
specimen,  with  a  clear  span,  L,  of  5.72  cm  (2.25  in.),  was  loaded  with  two  rollers  that  are  1.91  cm 
(0.75  in.)  apart  and  are  located  symmetrically  with  respect  to  the  center  of  the  specimen.  Loads  in 
the  samples  were  induced  by  turning  a  set  screw  located  along  the  centerline  of  the  strip.  The  mag¬ 
nitudes  of  the  loads  were  recorded  through  the  output  of  a  load  cell  located  at  the  tip  of  the  load 
screw.  The  maximum  stress  for  the  strip  is 

c,  =  3PbV(Lwt')  (19) 

where  P  is  the  total  load  measured  by  the  load  cell,  b  and  L  are  the  dimensions  shown  in  Figure  2, 
and  w  and  t  are  the  width  and  thickness,  respectively,  of  the  strip. 

Single-crystal  silicon  (Si)  strip  samples  normal  to  the  <00 1>  plane  were  fabricated.  The  sample 
dimensions  were  12.7  mm  (0.500  in.)  wide  and  0.394  mm  (0.0155  in.)  thick.  They  consist  of  two 
types  such  that  a  unidirectional  stress  could  be  applied  either  in  the  (1 10)  or  (100)  crystal  directions. 


4.2  Test  Apparatus 

In  the  Raman  stress  experimental  setup,  a  monochromatic  light  or  laser  shines  through  a  microscope 
onto  the  surface  of  a  specimen.  The  specimen  sits  on  a  translation  stage  whose  motion  is  controlled 
by  a  micrometer.  The  spectrometer  output  is  the  intensity  of  the  Raman  signal  as  a  function  of  fre¬ 
quency.  A  reference  spectrum  is  taken  for  an  unstressed  sample.  When  the  specimen  is  placed  in  a 
stressed  state,  the  spectmm  will  show  a  frequency  shift,  which  is  linear  to  the  applied  stress.  A  Ne 
penray  lamp  was  used  for  calibration  lines  to  reliably  and  accurately  measure  the  silicon  transition. 

The  measurements  were  conducted  using  an  argon  laser  line  (wavelength  488  nm).  The  amount  of 
frequency  shift  depends  upon  the  stress  state  and  can  be  determined  from  Eq.  (2).  The  diameter  of 


roller  roller 


Figure  2.  Schematics  of  strip  specimen. 
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the  laser  spot  is  typically  1  to  2  |j,m.  Therefore,  a  stress  variation  can  be  mapped  out  accurately. 
Raman  spectra  were  acquired  in  integration  mode  using  10  scans  at  2  s  per  scan. 


4.3  Results 

Figure  3  depicts  the  Raman  spectra  of  a  wafer  cut  from  an  annealed  single-silicon  crystal  boule.  The 
spectra  were  curve^fit  using  a  least-squares  fitting  procedure.  In  general,  the  fitting  error  was  on  the 
order  of  0.005  cm  and  is  too  small  to  be  considered  significant.  The  measured  average  center 
location  wave  number  is  520.28  cm  '  with  an  uncertainty  of  0.06  cm  '.  This  uncertainty  corresponds 
to  an  equal  biaxial  stress  of +14.7  MPa  (+2.1  ksi),  which  indicates  the  resolution  of  the  technique. 
Measurements  were  also  taken  on  an  annealed  polysilicon  layer  deposited  on  a  SiO^  layer  on  top  of  a 
single-crystal  silicon  substrate  (see  Figure  4).  These  measurements  showed  a  spatial  variability  of 
+38.6  MPa  (+5.6  ksi),  well  in  excess  of  the  measurement  resolution,  indicating  the  presence  of 
residual  stress.. 

A  calculation  was  made  based  on  the  p,  q,  r,  and  elastic  constants  of  single-crystal  silicon'.  For  sin¬ 
gle  crystal  silicon,  p  =-1.43  x  10“(s'%  q  =  -1.89  x  10“(s'^),  and  r  =  -0.59  x  10“(s'^).  The  calculated 
results  also  show  the  Raman  shifts  in  cm  '  as  follows: 

For  uniaxial  stress,  c:  =  7.83xl0"a  and  Af2j=  Af23=  2.00xl0"'a 

For  equal  biaxial  stress,  a:  AQ,=  AQ.^=  2.79xl0"o  and  4.00xl0'“’a, 

where  a  is  in  dyne/cm^ 

Preliminary  Raman  frequency  measurements  for  (100)  single-crystal  silicon  strip  with  (110)  applied 
unidirectional  stress  are  presented  in  Figure  5.  This  sample  and  its  stress  orientation  correspond  to 


Figures.  Raman  plot  for  unstressed  single-crystal  Si.  The  various  ciuves 

represent  single  spectra  sampled  at  different  locations  on  the  wafer. 
The  two  lower  amplitude  cases  were  obtained  without  refocusing. 
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Figure  4,  Raman  plot  for  polysilicon  deposited  on  Si02  on  single-crystal  Si. 
Again,  the  various  curves  are  spectra  obtained  at  different,  but 
randomly  located,  positions  on  the  wafer. 


Figure  5.  Frequency  shift  as  a  function  of  applied  stress  in  the  (1 10)  direction.  The 
sample  is  a  strip  of  (001)  single-crystal  silicon.  The  fitted  line  has  a  slope 
of -1.62  X  10’*^^  with  an  error  of  0.23  xlO'^^ 

the  equal  biaxial  case  discussed  above.  [However  the  stress  plotted  in  Figure  4  is  twice  that  for  the 
biaxial  case  above.]  The  measured  slope  is  -1.66  ±  0.23  x  10  ‘“(cn)/dyne)  compared  to  the  expected 
value  of  2.0  x  lO  '®.  The  agreement  is  within  20%  of  theory.  The  site  from  which  spectra  were 
measured  is  approximately  0.001  in.  in  size.  Spatial  variation  effects  (errors)  are  expected  to  be 
much  smaller  than  0.06  cm  '.  Additional  error  due  to  calibration  error  of  the  load  cell  is  estimated  to 
contribute  a  10%  error  in  the  stress  estimation.  The  Raman  frequency  measurements  with  (100) 
applied  unidirectional  stress  are  presented  in  Figure  6.  The  measured  slope  is  -7.95  x  10  ".  The 
theoretical  slope  is  -7.68  x  10  ".  The  agreement  is  well  within  experimental  errors.  Again,  measure- 


ments  were  performed  at  a  single  location  within  a  spot  of  0.001  in.  Deviation  from  the  fitted  line  is 
well  within  0.002  cm  ,  probably  representing  spatial  variations.  An  estimated  10%  error  in  stress  is 
due  to  a  calibration  error  on  the  load  cell.  Further  work  to  refine  the  measurements  and  reduce 
experimental  uncertainties  is  in  progress.  This  type  of  data  provides  the  necessary  calibration  and 
confirmation  of  the  theoretical  and  experimental  approach  taken  in  this  report. 
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5.  Conclusions  and  Discussion 


We  have  demonstrated  analytically  that  the  two-dimensional  surface  stresses  normal  to  the  <001> 
plane  of  a  cubic  crystal  have  a  unique  relationship  with  the  Raman  frequency  shifts.  In  other  words, 
once  the  Raman  frequency  shifts  of  a  cubic  crystal  (of  which  its  spring  constants  p,  q,  r  are  known) 
are  measured,  then  the  three  associated  surface-stress  components  can  be  uniquely  determined.  We 
have  also  demonstrated  that  the  Raman  frequency  shifts  measured  by  the  use  of  our  in-house  system 
under  uni-directional  stress  in  both  (1 10)  and  (100)  directions  agree  with  analytical  predictions. 

It  is  to  be  noted  that  values  of  p,  q,  and  r  for  single-crystal  silicon  were  used  in  all  the  calculations. 
The  values  for  polysilicon  may  vary  from  values  of  single-crystal  silicon  and  may  need  to  be  cali¬ 
brated.  In  fact,  because  polysilicon  tends  to  be  isotropic  in  microstructure,  there  may  only  be  two 
constants  needed  to  describe  the  Raman  behavior,  such  as  the  two  Lame  elastic  constants  for  an  iso¬ 
tropic  material.  An  experimental  exercise  is  required  to  address  this  issue.  Many  investigators  have 
used  the  p,  q,  and  r  values  of  single-crystal  silicon  for  polysilicon.  They  also  assumed  incorrectly 
that  the  stress  is  equal  biaxial  in  the  vicinity  of  discontinuities  such  as  structures.  For  the  backscat- 
tered  mode  used  in  many  of  these  studies,  it  is  not  possible  to  observe  the  effect  of  shear  stress  on  the 
Raman  spectra.  Shear  stress  is  expected  in  the  vicinity  of  long-aspect-ratio  structures  typically  used 
in  microelectronics. 

At  present,  our  Raman  system  uses  only  backscattered  light,  which  is  sensitive  to  only  one  of  the 
three  possible  frequency  shifts.  More  than  one  observation  geometry  and  polarization  condition  is 
required  to  detect  other  modes  or  solutions  to  Eq.  (2).'’  ’ 
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6.  Recommendations 


Since  Raman  spectroscopy  can  provide  very  localized  stress  measurements  for  mini  to  micro-sys¬ 
tems  such  as  microelectronics  and  MEMS  systems,  it  has  a  strong  potential  as  an  aid  in  the  estab¬ 
lishment  of  reliability  assessment  techniques.  However,  although  its  application  to  single-crystal 
systems  is  well  understood,  there  is  a  need  to  expand  our  understanding  to  polycrystal  systems.  It  is, 
therefore,  recommended  that  the  Raman  technique  in  stress  assessment  for  polycrystal  systems  such 
as  polysilicon  be  investigated. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  specializing 
in  advanced  military  space  systems.  The  Corporation’s  Laboratory  Operations  supports  the  effective  and 
timely  development  and  operation  of  national  security  systems  through  scientific  research  and  the  application 
of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staff’s  wide-ranging  expertise 
and  its  ability  to  stay  abreast  of  new  technological  developments  and  program  support  issues  associated  with 
rapidly  evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and  CCD 
detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid  state  laser 
design,  micro-optics,  optical  communications,  and  fiber  optic  sensors;  atomic  frequency  stan¬ 
dards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam  control, 
LIDARyLADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery  electro¬ 
chemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and  process¬ 
ing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites;  development  of 
advanced  deposition  processes;  nondestructive  evaluation,  component  failure  analysis  and  reli¬ 
ability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion;  analysis  and  evaluation 
of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle  fluid  mechanics,  heat  trans¬ 
fer  and  flight  dynamics;  aerothermodynamics;  chemical  and  electric  propulsion;  environmental 
chemistry;  combustion  processes;  space  environment  effects  on  materials,  hardening  and  vul¬ 
nerability  assessment;  contamination,  thermal  and  structural  control;  lubrication  and  surface 
phenomena. 

Space  Science  Application  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric  physics, 
density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmospheric  radiation; 
solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  surveillance,  imaging, 
remote  sensing,  and  hyperspectral  imaging;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and  magnetosphere;  effects  of  elec¬ 
tromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation,  design  fabrica¬ 
tion  and  test;  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions,  atmos¬ 
pheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of  mis¬ 
sile  plumes. 

Center  for  Microtechnology:  Microelectromechanical  systems  (MEMS)  for  space  applica¬ 
tions;  assessment  of  microtechnology  space  applications;  laser  micromachining;  laser-surface 
physical  and  chemical  interactions;  micropropulsion;  micro-  and  nanosatellite  mission 
analysis;  intelligent  microinstruments  for  monitoring  space  and  launch  system  environments. 

OlBfice  of  Spectral  Applications:  Multispectral  and  hyperspectral  sensor  development;  data 
analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral  imagery  to 
defense,  civil  space,  commercial,  and  environmental  missions. 


